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INTRODUCTION 


The use of acoustic liners to reduce aircraft engine noise has been the 
subject of intense investigation by both governmental and industrial groups for 
over a decade. Figure 1 is an illustration of three types of acoustic liners 
which have potential noise suppression benefits. Earlier research was concen- 
trated on uniform liners. Recent efforts have centered on axially segmented 
liners where the liner admittance changes with distance along the duct. There 
has been little previous work on peripherally segmented liners. 

Certain characteristics of axially segmented liners may preclude their 
application to aircraft engine noise suppression. The optimum design of axially 
segmented acoustic liners requires an accurate description of the modal distri- 
bution of acoustic energy; however, the number of modes generated in a typical 
aircraft engine is large, and the problem of determining the modal amplitude and 
relative phasing of these source modes is difficult (if not impossible) . Engine 
nacelles of today * s aircraft may not provide sufficient length for the effective 
use of these liners. Further, axially segmented liners perform best in the 
lower range of the engine noise spectrum. They do not perform as well in the 
higher frequency range where greater suppression is needed. 

The peripherally segmented liner has a wall admittance which changes in the 
circumferential direction. These liners may have significant advantages when 
compared to uniform or axially segmented liners. The duct length is less impor- 
tant because of the fact that segmentation is around the duct circumference 
instead of along the axis. This is important since high-bypass-ratio engines 
have inlet ducts with diameters as large as their lengths. Peripherally seg- 
mented liners redistribute acoustic energy of a circumferential wave into 
higher-order circumferential waves and, since these higher order waves are more 
rapidly attenuated, an improved suppression may result. 

Mani (ref. 1) tested a peripherally segmented liner and found that some 
additional suppression beyond that achieved with a uniform liner could be 
observed over a broad frequency range. On the basis of this work, he was 
granted a patent on the concept. However, there has been no extensive experi- 
mental testing of this idea. In addition, no theoretical analyses of such 
liners are available. The separation-of-variables technique cannot be used to 
obtain the solution for the acoustic field in ducts with peripheral liners. 

Thus, the solution for the acoustic field in the duct's cross section is more 
complex than for uniform or axially segmented liners, where separation of vari- 
ables is possible in the circumferential direction. 

Watson (ref. 2) attempted a multimodal, finite-element solution to this 
problem. However, the number of peripheral strips was severely restricted 
because of matrix size limitations, and the effects of a realistic sound source 
could not be evaluated. Abrahamson (ref. 3) attempted a direct numerical solu- 
tion of the linearized three-dimensional Euler equation with variable-admittance 
boundary conditions. He concluded that this method could not succeed using 
present-day computers because of the size and character of the coefficient 
matrix. Astley et al . (ref. 4) also attacked the problem using a multimodal. 



finite-element method, including the effects of an arbitrary mean flow. 

Astley's approach is sxibject to the same limitations as the method of refer- 
ence 2- Thus, there is no available method for computing the acoustic field in 
ducts with peripherally segmented liners. Such a method is necessary before an 
evaluation can be made of the effectiveness of circumferential liners as sound 
attenuators. A method is developed herein which can account for the sound 
source and a large number of peripheral strips. The primary objectives of this 
paper are to compute the eigenvalues and eigenfunctions in ducts with peripher- 
ally segmented liners and to compare their properties to those of the uniform 
liner. It is shown how redistribution of the acoustic energy of a single cir- 
cumferential wave into a multitude of other circumferential waves is a basic 
physical property of peripherally segmented liners which is not possible for 
uniform or axially segmented liners. A secondary objective of this paper is to 
demonstrate the use of uniform-liner eigenfunctions to calculate the eigen- 
values and eigenfunctions of peripherally segmented liners. Finally, attenua- 
tion characteristics of a peripheral liner are compared to those of an optimal 
uniform liner for plane-wave sources. 


SYMBOLS 

a duct radius 

A amplitude of pressure wave 

B admittance matrix 

c speed of sound 

i unit imaginary mamber 

( ) imaginary part of complex expression 

Bessel function of first kind with order m 
K wave number 

L predetermined length over which transmission loss is computed 


M 

number 

of 

m values. 

m = 0, 

1, 

1 — 1 

N 

number 

of 

n values. 

n = 0, 

1, 

, . . , N-1 


normalization constant for basis function 
P(R,0,Z) acoustic pressure field 

P£(R,0) acoustic pressure eigenfunctions 

P£j coefficients of basis functions in series acoustic pressure 

eigenfunctions 
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Q(R/0) source pressure 

Qj coefficients of basis functions in series for source pressure 

r, 7, dimensional coordinates 

R,Q,Z cylindrical coordinates 

Re ( ) real part of complex expression 

S number of peripheral strips 

[s/N] integer division of integer s by integer N 

T periodicity of peripherally segmented liner 

3(0) acoustic admittance of duct wall 

6j Kronecker delta function 

X eigenvalue 

p ambient density of medium 

(t)j (0) basis functions for series expansion of 3(9) 

lpj(R,0) basis functions for series expansions of P^{R,Q) and Q(R,0) 

Laplace operator 
Subscripts : 

j basis function index 

£ eigenfunction index 

m circumferential harmonic number 

n radial (basis) eigenfunction index 

s, t integers 

A bar over a symbol denotes a basis function or a constant- A prime 
denotes the first derivative with respect to distance. 

STATEMENT OF PROBLEM 

Consider a semi-infinite circular duct with radius a, as shown in fig- 
ure 2. A lining material with specific acoustic admittance 3(9) is placed 
along the inner wall, and an initial noise source (such as one generated by an 
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aircraft engine) is given at z = 0. The purpose of the acoustic liner is to 
attenuate the noise as it propagates along the duct. The mathematical problem 
is to determine the acoustic field inside the duct. The effectiveness of the 
lining material may be evaluated from the solution for this field. 

Acoustic waves are governed by the Helmholtz equation 


V^P(R,0,Z) + K^P(R,0,Z) = 


( 1 ) 


in which a time dependence of the form e has been assumed, K = is the 

c 

dimensionless wave number, P is the dimensionless acoustic pressure, a) is 
the angular frequency, t is time, and c denotes the ambient speed of sound. 
In equation (1) , all distances are referred to the duct radius a and the 
acoustic pressure is referred to the quantity pc^ in which p is the ambient 
density of the medium. 

The acoustic boundary condition along the outer wall (R = 1) is expressed 
in terms of the admittance of the lining material 3(0) f 


dP (R,9,Z) 
dR 


= iK 30) P(R,0,Z) 


(2) 


At the entrance plane, the boundary condition is assumed known as 


P(R,0,O) = Q(R,0) 


(3) 


Capital letters are used for dimensionless groups throughout this paper. 

Equations (1) to (3) constitute a boundary value problem for the acoustic 
pressure field P(R,0,Z). A solution must be obtained before the effectiveness 
of the peripheral liner can be evaluated. Analytical solutions to these equa- 
tions are not available since the acoustic admittance 3 depends on the azi- 
muthal coordinate 0. Previous attempts at finite element solutions have had 
limited success because of the size and character of matrices used in the 
analysis. Thus, an approximate series expansion is used in this paper. 


EIGENFUNCTION EXPANSION TECHNIQUE 

The solution to equations (1) to (3) is expressed in the form (ref. 2) 


P(R,0,2) 


2MN-N-1 


2 

£=0 


Pj^(R,0) e 


iK^Z 


(4) 
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in which each eigenfunction P^CR^B) satisfies a Helmholtz equation 


y2 Pj^(R,0) 


+ xl Pjj^(R,0) = 0 


(5) 


with the homogeneous boundary conditions 


9Pj^(r,0) 

9r 


+ 


iK 8(6) P£(R,6) 


0 at R = 1 


( 6 ) 


where 

value 


is the two-dimensional Laplace operator in R and 0. The eigen- 
is 


A^ - 


(7) 


The eigenfunctions satisfying equations (5) to 
This can be proved using Green's theorem 


(7) are also orthogonal. 



V^P, 


m 


m 




dA 




P„ 9p 

£ m 

0R 


P 3p.\ 
m £ I 
9 r / 


dc 


( 8 ) 


In view of the fact that both P^^ and P^ satisfy equations (5) and (6) , 
equation (8) becomes 


(4 - ^m) = 0 


(9) 


If ^ equation (9) becomes 


I 




dA = 0 


(10) 


where m ^ £. Further, the eigenfunctions will be required to form an ortho- 
normal set so that 



( 11 ) 


Equations (5) to (7) constitute an eigenvalue problem which must be solved 
before the acoustic field of equation (4) can be determined. In this paper, a 
numerical solution to this equation is obtained using a second modal expansion 
for the eigenfunctions Pj^(R,0). 
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SERIES SOLUTION FOR EIGENFUNCTIONS 


To solve the eigenvalue problem for a peripherally segmented liner, a 
Galerkin method is employed (see ref. 5) which simultaneously considers the 
field equation (eq. (5)) and the admittance boundary condition (eq. (6)). 
Solutions of the form 


2MN-N-1 

P£(R,0) == 

s=0 

are sought in which 


IPg (R,0) 



COS m0 


sin m0 


[s/n] 

(S - N - MN)/n] 


n = 



for s < MN 
for s - MN 

for s < MN 
for s - MN 

for s < MN 
for s - MN 


( 12 ) 


(13) 


(14) 


(15) 


in which the notation [s/n] denotes integer division of the integer s by 
the integer N. The normalization constant Njj^j^ is determined so that the 
basis functions i|^g(R,0) are orthonormal 


Nmn 

mn 


= jj Jm( Vn^) “9) 


dA 


(16) 


and the eigenvalues are uniform-liner eigenvalues chosen to satisfy the 

transcendental equation 



iK3 

m\ mn/ 


0 


(17) 


6 



The constant admittance B may be arbitrarily chosen. Two of the possible 

1 

choices are 3 = 0, for a hard-wall duct, and 3 = \ 3(9) <3.0, the average 

value of 3 for the peripherally segmented liner. The use of these two basis 
values will be discussed in the following sections. 

The unknown coefficients are determined by Galerkin's method. 

Galerkin’s condition is 






where t = 0, 1, ..., 2MN-N-1. Integrating equation (18) by parts to transfer 
the Laplace operator to the basis functions gives 


jy (^\ + = y 




Eliminating V and 3Pn/9R from equation (19) gives 


(4 - dA = iK J 


(3 - 3)P£4^t dc 


and substituting the expansion for into equation (20) gives 

2MN-N-1 


l^£t = 5! (‘^st^t ®ts)Pji 


where 




(3 - 3)'Po’I^ 4- dc 


Equation (21) determines the vector P^^ to within an arbitrary constant. 
This constant is defined through the normalization condition 


2MN-N-1 


1 


Equations (21) to (23) represent a system of linear equations which can be 
solved for the eigenvalue and series coefficients P^^ using standard 

matrix techniques. In this paper, the QR algorithm is employed to obtain the 
solution to this equation (ref. 6) . 
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CHOICE OF PERIPHERAL LINER AND SOURCE 


The source pressure function is usually generated in a hard-wall section. 
Thus, the source will be specified in terms of hard-wall duct modes. The source 
pressure function Q(R,0) is expanded in the form 

2MN-N-1 

Q(R,0) = ^ 2s IPs 

s=0 


in which the hard-wall basis functions are obtained from equation (17) by 

setting 3=0. 

One may also assume that the fabrication of a peripheral liner 3(9) 
will consist of some sort of periodicity as shown in figure 3. It is assumed 
that the admittance of the liner repeats after each A0 on the interval 

2tt 

0 ^ 0 ^ 2TT where A0 = — . A liner of this type is referred to as a liner 

with periodicity T, where T denotes the number of times the liner repeats on 
the interval 0 ^ 0 ^ 277 . Within each interval, the liner may be subdivided 
into strips which do not form a repeating pattern. 

Because of its periodicity, the admittance function 3(9) is expanded into 
the fourier series 


3(0) = 2 ^25) 

s=0 

where 


^ cos mT0 

for s even 




(26) 

Vsin mT0 

for s odd 


m = [ (s + l)/2] 


(27) 


EFFECTIVENESS OF LINING MATERIAL 


Before proceeding with the solution to equations (5) to (7) , it is useful 
to develop an expression for the attenuation produced by the liner in terms of 
parameters which have been introduced. The axial acoustic intensity at any 
axial position in the duct is 


I(R/0/Z) 



(28) 
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in which the nondimens ional acoustic intensity is referred to pc , Re [ ] 
denotes the real part of the complex expression enclosed within the brackets, 
and the superscript asterisk indicates the complex conjugate. The total 
acoustic power is the integral of the acoustic intensity across the cross 
section 


W(Z) 


II 


I dA 


( 29 ) 


and the transmission loss, or decrease in level of the acoustic power from 
Z = 0 to Z = L, can be written as 


TL = 10 logj^Q [w(0)/(W(L^ (30) 

Generally, the specific admittance 3 is chosen so as to maximize the trans- 
mission loss. 

In this paper, the eigenmode p£ is expanded in terms of the hard-wall 
basis functions. These basis functions were chosen not only because the eigen- 
values are real and appear in standard tables but because the resulting 

expression for the acoustic power reduces to 


2MN-N-1 2MN-N-1 
W(Z) . ^ j 

^=0 ^=0 

2MN-N-1 


"HZ 




2 

s=C 

-N- 

2 


^Zs^Zs 


s=0 

2MN-N-1 


: a' 


K Re 1 ApA'p Kp I n 'p e 

A/ A/ A/ A/ A/ 


i(Kj-K^)Z 


(31) 


(32) 


(33) 


s=0 


The function W(Z) cannot be simplified to a double summation if 3^0. 


RESULTS AND DISCUSSION 

In this section, eigenvalues and eigenfunctions are presented for uniform 
and peripherally segmented liners. Emphasis is placed on the solution for the 
eigenvalue and the eigenfunction P£(R,0) for peripherally segmented 

liners. The way these solutions change with the number of peripheral strips and 
the wave number K of the system is of fundamental interest. First, confidence 
in the analysis is obtained by demonstrating that it agrees with exact results 
which are possible when the liner admittance is independent of the coordinate 0. 
Afterward, results are presented for peripherally segmented liners. 
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Uniform Liners 


Analysis When the admittance 6 is independent of the coordinate 0, an 
exact solution to equations (5) to (7) can be obtained by requiring that 
3=3- Equation (21) becomes 



where s = 0, 1, 2MN-N-1. The exact solution to this equation, with 

the normalization condition, is 



where m and n are related to s by equations (14) and (15) . Observe that 
the solution Pjj^(R,0) has the following properties for the uniform liner: 

(1) Each eigenfunction is either an even or odd function of 0 

(2) The eigenvalues for the even and odd eigenfunctions are identical 

(3) Each eigenfunction is composed of only a single function of m0 

These observations are important, since it will be shown in the following sec- 
tion that the eigenfunctions for a peripheral liner will have_none of these 
properties in general. Further, if 3 is chosen such that 3 3f then, 

^0 in general and equation (21) cannot be solved exactly (the solution 
was obtained employing the QR algorithm) . However, the solution for the eigen- 
functions P^(R,0) and eigenvalue will have the three properties dis- 

cussed earlier. These three properties are_preserved because the matrix equa- 
tion which governs equation (22) when 3 3 is block diagonal, with each 

diagonal block giving the solution for either an even or odd eigenfunction. 

This eigenfunction is composed of only a single function of m9 . 

Now consider the eigenfunction expansion of the acoustic field as given by 
equation (4) . The modal coefficient is 

pp 2M^N-l 

= U Q(R,0) P^(R,0) dA = 2 

s=0 

Now for uniform liners, A^ is zero for any standing wave (cos m0 or sin m0) 
not present at the source. Thus, for uniform liners, the acoustic energy in an 
even standing wave of the source remains in that even standing wave in the lined 
section. This is also true for odd standing waves. It is not possible, there- 
fore, to redistribute the energy of a single standing wave of the source into a 
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multitude of standing waves in the lined section with a uniform liner. In 
addition the uniform liner restricts the energy in even standing waves of the 
source to even waves in the lined section. Similarly, the energy in odd stand- 
ing waves of the source is restricted to odd standing waves in the lined sec- 
tion. In the next section, it is shown that peripherally segmented liners do 
allow the redistribution of acoustic energy out of even standing waves of the 
source into even and odd standing waves in the lined section. It is also shown 
that peripherally segmented liners allow the redistribution of acoustic energy 
out of a given standing wave of the source into a multitude of standing waves 
in the lined section. These properties of peripherally segmented liners may 
make them better attenuators of sound than uniform liners. 

Ei genvalue computation .- Attention is now focused on the eigenvalues com- 
puted from the present analysis and how they approximate the exact eigenvalues 
satisfying equation (17). In order to compare the two solutions, equation (17) 
has been solved using a Newton-Raphson iterative procedure (ref. 7) with 3 = B. 
Results for uniform liners will be restricted to optimum liners for plane-wave 
sources where the optimum values of B have been taken directly from the 
results of Lester and Posey (ref. 8) . All uniform liner results are for 
m = 0, 1, 2 and N = 10. 

Table I allows comparison of eigenvalues computed from equation (21) to 
those computed from solving the exact transcendental equation for m = 0 with 
K = 2- Note that the present analysis computes the eigenvalues with the larger 
magnitudes more accurately than those with the smaller magnitudes. This is not 
surprising since lower-order eigenvalues can be expected to have more simple 
mode shapes than higher-order ones. Thus, it can be expected that more terms 
are needed in a Bessel function series to represent functions which are closer 
to plane waves. Overall, the agreement between the exact eigenvalues and those 
of the present analysis ranges from one to four significant figures. 

Eigenvalues are presented in tables II and III for m = 1 and m = 2, 
respectively. Values of K and B a.re identical to those of table I. The 
trends of these tables are consistent with those of table I . 

Results are presented in table IV for a nondimensional wave number of 
K = 5 for m = 0. Eigenvalues computed from equation (21) compare favorably 
with those computed from the exact transcendental equation. Computations for 
m = 1 and m = 2 are given in tables V and VI, respectively. The value of 
the admittance B is identical to that of table IV. These results show that 
better comparison is obtained at this higher value of K. The eigenvalues com- 
puted in tables IV, V, and VI agree to within two to three significant figures 
with those of the exact transcendental equation. 

Computations have also been obtained for K = 10 . These results are shown 
in tables VII, VIII, and IX for m = 0, 1, and 2, respectively. Results in 
these tables are consistent with those of the previous tables. The highest 
order eigenvalues are computed more accurately than the lower order ones. 

Again, agreement between the eigenvalues computed from the exact equation and 
those of the present analysis is within two to three significant figures. 
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Overall, the agreement between the eigenvalues tends to improve with increasing 
wave number K and mode order £• 

Eigenfunctions > - Attention is now focused on the characteristic functions 
P^(R,Q) computed from the analysis presented in this work. As in any accept- 
able approximate formulation, the series solution computed here must converge 
in some sense to the exact solution of the problem. Because basis functions 
used here are for 3=0, the expansion given by equation (12) cannot possibly 
satisfy the boundary conditions at the outer wall. It is useful to study the 
effects of this choice of basis functions on the approximate solution. Further, 
since the circiimf erential dependence of each eigenfunction will match for a 
uniform liner, only the radial dependence of the eigenfunctions will be studied. 

The real and imaginary parts of the £ = 0 mode for K = 2 are plotted 
in figures 4 and 5, respectively. Agreement between this analysis and the 
exact analysis is poor for this lowest order mode- Results for the £ = 4 and 
£ = 9 modes are shown in figures 6 to 9 for K = 2. The results for these 
higher-order modes are virtually indistinguishable from the exact results, even 
at the boundary- In general, the higher the mode number, the more accurately 
this analysis computes the eigenfunction. 

The real and imaginary parts of the £ = 0 mode for K = 10 are plotted 
in figures 10 and 11, respectively. At this higher value of K, this lowest- 
order mode is in good agreement with the exact value. Results for the £ = 4 
and £ = 9 modes are given in figures 12 to 15 for this higher value of K. 
Results computed from this analysis are identical to the exact values for these 
higher-order modes. 

Characteristic functions have also been computed for K=3, 4, 5, 6, 7, 
and 8 for various values of £. Results for these will not be presented for 
the sake of brevity. However, trends show that the accuracy of the modes com- 
puted from this analysis increases with increasing K, with the lowest-order 
mode being less accurate than higher-order ones. 

Attention is now focused on the radial derivative of the characteristic 
function P£(R,0). How the radial derivative computed from this analysis com- 
pares with the exact value is of special interest because of the fact that the 
hard-wall duct modes are used as basis functions. Recall that the radial 
derivatives computed from this analysis cannot compare with the exact value at 
the boundary unless 3=0, or P£(1,0)=O. 

The real and imaginary parts of Sp^/Sr for the £ = 0 mode are plotted 
in figures 16 and 17, respectively, for K = 2. Overall, values computed from 
this analysis do not compare to the exact values for this lowest-order mode. 
Further, results computed from this analysis are less accurate as the outer 
boundary (R = 1) is approached. 

Results for the real and imaginary parts of 3P£/9 r for the £ = 4 and 
£ = 9 modes are plotted in figures 18 to 21 for K = 2. Results for these 
higher-order modes are virtually indistinguishable from the exact values, with 
the exception of the outer boundary (R = 1) . 
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Shown in figures 22 and 23 are results for for £ = 0 with 

K = 10. Observe that the analysis used here computes the radial derivative 
of P£ more accurately at this higher value of K than at K = 2. Overall, 
this lower^order mode has not been accurately computed by this analysis at this 
higher K, and the accuracy decreases as the boundary R = 1 is approached. 

Results for the real and imaginary parts of Bp^/Sr for the £ = 4 and 
£ = 9 mode are plotted in figures 24 to 27 for K = 10 - Results for these 
higher-order modes are in good agreement with the exact values except near the 
outer boundary. 

Transmis sion losses .- In this section, the transmission losses predicted 
from equations (30) to (33) for optimum uniform liners are compared to results 
predicted from the analysis of reference 8. The source pressure distribution 
was assumed to be a plane wave so that O(R,0) = 1. The duct geometry was 
chosen so that the length to diameter ratio was unity (L = 2) . 

Transmission losses are presented in figure 28, The transmission loss 
spectrum is given for a frequency range 2 - K - 20. The losses predicted from 
the analysis of this work are identical to those of the method of reference 8 
for K > 4. However, for K < 4, the transmission losses predicted from this 
analysis are less than those predicted in reference 8. As the frequency 
decreases, the difference between the two curves in figure 28 increases. 

The discrepancies between the transmission losses at low frequencies 
(K < 4) predicted from this analysis and from that of reference 8 are not sur- 
prising. Thus, for K < 4, only the £ = 0 mode is cut on in the duct 
(ref. 8) . However, the analysis of this work cannot accurately resolve the 
£ = 0 mode at low frequencies . 


Peripherally Segmented Liners 

Analysis .- Attention is now focused on peripherally segmented liners. To 
understand why the eigenfunction P£(R, 0 ) for a peripheral liner has properties 
fundamentally different from those of a uniform liner, the basic structure of 
equation (21) is studied when 3 is a function of the coordinate 0. In this 
instance, the matrix of coefficients 6^3 is not block diagonal as for the uni- 
form liner. As a consequence, the energy in either an even or odd mode of the 
source is redistributed into other even and odd modes in the lined section. 

This also results in the generation of a multitude of standing waves for each 
eigenfunction P£(R,0). As a consequence, the energy in a single standing wave 
of the source is redistributed into a multitude of standing waves in the lined 
section of the duct. The results of this section will show that this does 
indeed occur. 

Now the number of unknowns in equation (21) will be large. For example, 
the order of this matrix equation will be 190 if M = 10 and N = 10. However, 
if the function 3 ( 0 ) is restricted to an even function, = 0 for 

s > MN and t ^ MN. As a result, the eigenfunction P^^(R,0) will be composed 
of either even or odd modes. Further, if the source produces only even modes. 
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only the even functions need be included in equation (12) and the order of the 
system is reduced by (MN - N) . Thus, results in this section are restricted to 
admittance functions which are even and, in addition, the source is assumed to 
be composed of only even modes. Results of the previous section also show that 
10 radial modes (N = 10) appear to be sufficient in the series expansion of 
equation (12) for uniform liners. The number of radial modes has been termi- 
nated at N = 10 for peripherally segmented liner results. Note also that up 
to this point there is no rule for choosing the number of standing waves in the 
expansion given by equation (12). Finally, since 3(9) is assumed an even 
function, 3 (0) = 3 (277 - 0) and the form of 3(0) need only be specified on 
the interval (0,77) . Results for peripheral liners are restricted to a liner 
having a periodicity of one (T = 1) in which the admittance segments alternate 
between 3i and 32 with a given number of strips S on the interval (0,77) , 
as shown in figure 29. 

Eigenvalues Table X contains eigenvalues for a two-strip liner (S = 2) 
for various values of M. The liner admittance has been taken as 
3]_ = 1.54 - 1.29i and ^2 ~ 0*9 O-Oi with the dimensionless frequency K 
chosen as 2. These results show that the eigenvalues are significantly differ- 
ent for the lower values of M, but do converge as M increases. In fact, 
there is very little change in the first 10 eigenvalues when M is increased 
from 6 to 10. Table XI contains eigenvalues for K = 5, with all other param- 
eters identical to those of table X. The effects of M at this higher value 
of K are consistent with those at the lower value of K. Eigenvalues have 
also been computed for K = 4, 6, 8, and 10, and the effects of increasing M 
were always consistent with those of tables X and XI. 

The effects of varying the number of strips on the eigenvalues of periph- 
eral liners for K = 2 and K = 5 are shown in tables XII and XIII, respec- 
tively. Values of 3]_ and 32 employed in the two tables are identical to 
those of table X, and 10 standing waves have been employed (M = 10) . It is 
clear from these two tables that the number of strips affects the eigenvalues 
and subsequently the eigenmodes P£(R,0) even though the amount of the duct 
which is acoustically treated for each configuration is still 50 percent. 

Thus, the attenuation properties of peripherally segmented liners will depend 
on the number of peripheral strips, ^yen if the amount of lined section is kept 
fixed. 

Eigenfunctions . - The analysis presented in this work can also be used to 
investigate the eigenfunctions of peripherally segmented liners. Valuable 
insight about the eigenfunctions can be obtained by studying the coefficients 
P£s which go into the series expansion for the eigenfunction P£(R,0) in 
equation (12) . In table XIV the coefficients P^^ computed for the £ = 0 
mode of table XII are shown. Only values of £ which correspond to n = 0 are 
shown (i.e., each value of s corresponds to a different order standing wave). 
The first column of coefficients in this table is for the uniform liner for 
which 3-j_ = 32 = 1-54 - 1.29i and for the same value of K as the peripheral 
liner. Results of table XIV confirm that the peripheral liner generates a 
multitude of standing waves regardless of the number of segments. This is 
borne out by the fact that P^^ 0 for any s for the peripheral liners . 

Contrast this to the uniform liner coefficients in which Pn = 0 for s 0. 
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Results similar to those of table XIV have been computed for other values of £ 
in addition to £ = o and for other values of K, and ^or each 

value of £, the peripheral liner eigenfunction P£(R, 0 ) consisted of a multi- 
tude of standing waves, while the uniform liner contains only one standing 
wave. Thus, the multimodal expansion of the acoustic field for peripherally 
segmented liners (eq. ( 4 )), will contain a multitude of standing waves for any 
given source . 

Tr ans mis sion losses . - Baseline trends for the transmission loss of a 
peripherally segmented liner are given by the attenuation of the least atten- 
uated mode (this would correspond to setting = 0 in eq. ( 33 ) for £ 7 ^ 0 ) . 

The attenuation of the least attenuated mode is -8.69 Im(KQL) dB in which 
Im ( ) denotes the imaginary part of the complex expression within the paren- 
theses and dB denotes decibels of sound. 

The attenuation of the least-attenuated mode £ = 0 is assumed to provide 
a rough estimate of the transmission loss for the peripheral liner. This 
attenuation would equal the transmission loss if the duct was infinitely long 
and the mode £ = 0 carried all the acoustic energy. Figure 30 allows com- 
parison of these attenuations for uniform and peripherally segmented liners 
for L = 1. The uniform liner is chosen so that the attenuation is optimum at 
K = 3 (this implies an admittance of 6 = 0.81 - 0.78i); whereas, the periph- 
erally segmented liner is chosen so that Bj|_ = 0.81 - 0.78i and 82 = 0- 
Results of figure 30 show that peripheral liners with a small number of strips 
(S = 2 and S = 4) have poor attenuations compared to the uniform liner. The 
S = 6 and S == 10 liners attenuate more for frequencies >3.5 than the uniform 
liner. This is rather surprising, since only 50 percent of the duct wall is 
treated with peripheral strips; whereas, the unifomn liner covers the entire 
wall. Peripheral liners with a large number of strips have broader attenuation 
characteristics than uniform liners, and the peak attenuation occurs at a higher 
frequency. These characteristics may be a significant advantage for actual 
engine sources which produce noise over a wide frequency range. 

Turning now to the effects of the sound source, figure 31 is a comparison 
of transmission losses for an optimum uniform liner to those of an arbitrarily 
chosen peripherally segmented liner. Results in this figure are for a plane- 
wave source (Q = 1) with L = 2. The uniform liner is optimum at K = 5 so 
that 3 = 0.34 - 0.59i for this liner (ref. 8 ). The admittances of the 
peripheral liner are chosen such that - 0.59i and S >2 ~ 0.5 - 0.2i. 

The transmission loss was computed from equations (30) to (33) . Both peripher- 
ally segmented liners give transmission losses comparable to that of the uniform 
liner for K S 6 and for K ^ 17. For values of K between these values, the 
peripheral liner gives greater transmission losses than the uniform liner for 
this plane-wave source. Further, the S = 4 peripheral liner gives a 
50-percent increase in the attenuation over the uniform liner for K = 8 . 


CONCLUSIONS 

An analysis method has been developed for computing the acoustic field and 
transmission losses in ducts with peripherally segmented liners, arbitrary 
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sources, and any number of peripheral strips. Based upon this work, the 
following conclusions are made: 

(1) Peripherally segmented liners redistribute the energy in standing 
waves composed of only a single circumferential mode at the source into other 
waves which contain a multitude of circumferential modes in the lined section. 

(2) The accuracy of the eigenvalues and eigenfunctions computed from this 
analysis increases with both the frequency and the mode order. 

(3) Transmission losses computed from this analysis were determined to be 
accurate for wave numbers greater than four, which is the approximate cut-on 
value of the first order radial mode in a hard-wall duct. 

(4) Results based on the attenuation of the least-attenuated mode show 
that peripherally segmented liners can attenuate more sound over a broader fre- 
quency range than a uniform liner, although 50 percent of the peripheral liner 
was hard wall . 

(5) Finally, results show that, for a plane-wave source, a nonoptimal 
peripherally segmented liner may attenuate as much sound as an optimized uni- 
form liner at the optimal point, while giving more noise suppression at other 
frequencies. These results lend credence to the concept of a peripherally 
segmented liner and encourage further analysis utilizing optimized peripherally 
segmented liners and more realistic sources. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 5, 1981 
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TABLE I.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 


(K = 2, 3 = 1.54 - l,29i) 


Radial 

number 

mode 

a 

Exact solution to 
equation: 

transcendental 
- 213 ° 

Numerical solution 
computed from solving 
equation (21) 

0 


2.885 - 

0.489i 

2.950 - 0 

.461i 

1 


3.062 - 

3.083i 

2.740 - 3 

.120i 

2 


6.592 “ 

.422i 

6.600 - 

.373i 

3 


9.902 - 

.299i 

9.903 - 

.265i 

4 


13.122 - 

.230i 

13.121 - 

.204i 

5 


16.310 ~ 

.186i 

16.309 - 

.165i 

6 


19.482 - 

.157i 

19.481 - 

.138i 

7 


22.645 - 

.135i 

22.644 - 

.117i 

8 


25.803 - 

.119i 

25.802 - 

.lOli 

9 


28.957 - 

.106i 

28.956 - 

.087i 


TABLE II.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE -WAVE SOURCE 

(K = 2, 3 = 1.54 - 1.291) 


Radial 

number. 

mode 

i 

Exact solution to 
equation: 

transcendental 

- 2i3 = 0 

Numerical solution 
computed from solving 
equation (21) 

10 


3.146 - 

3.0411 

2.860 - 3 

.0701 

11 


4.683 - 

.5341 

4.720 “ 

.4761 

12 


8.199 - 

.3591 

8.202 - 

.3201 

13 


11.473 - 

.2631 

11.472 - 

.2351 

14 


14.684 - 

.2071 

14.683 - 

.1851 

15 


17.869 - 

.1711 

17.868 - 

.1521 

16 


21.040 - 

.1461 

21.039 - 

.1291 

17 


’ 24.204 - 

.1271 

24.203 - 

.1111 

18 


27.362 - 

.1121 

27.361 - 

.0961 

19 


30.517 - 

.lOli 

30.516 - 

.0831 
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TABLE III.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 


(K = 2, e = 1.54 - 1.29i) 


Radial 
mimber , 

mode 

1 

Exact 

equation 

solution to 

transcendental 

- 213 "" ° 

Numerical solution 
computed from solving 
equation (21) 

20 



3-537 

- 

2.876i 

3.277 - 2 

.8601 

21 



6.203 


.4851 

6.216 - 

.4331 

22 



9.680 

- 

.3181 

9-680 -- 

.2861 

23 



12.962 

- 

.2381 

12-961 - 

.2141 

24 



16.183 

- 

.1911 

16.182 - 

.1711 

25 



19.377 

- 

.1591 

19.376 - 

.1421 

26 



22.555 

- 

.1371 

22.554 - 

.1211 

27 



25.725 

- 

.1201 

25.723 - 

.1051 

28 



28.888 

- 

.1071 

28.886 - 

.0921 

29 



32.046 

- 

.0961 

32.045 - 

.0921 


TABLE IV.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 

(K = 5, 3 = 0.34 - 0.591) 


Radial mode 
number , Z 

Exact solution to transcendental 
equation: A£ ~ ^ 

Numerical solution 
computed from solving 
equation (21) 

0 

3.058 - 0.2771 

3.096 - 0.2601 

1 

1.668 - 3. 4851 

1.476 - 3.3701 

2 

6.589 - .2171 

6.609 - .1941 

3 

9.881 " .1591 

9.894 - .1411 

4 

13.101 - .1241 

13.111 - .1091 

5 

16.291 - .1011 

16.299 - .0891 

6 

19.465 - .0861 

19.422 - .0741 

7 

22.630 - .0741 

22.637 - .0641 

8 

25.790 - .0651 

25.796 - .0551 

9 ! 

28.945 - .0581 

28.953 - .0471 
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TABLE V.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 


(K = 5, e = 0.34 - 0.591) 


Radial 

nimber. 

mode 

£ 

Exact solution to 
equation: 

transcendental 
- 51$ Ji(Aj^) = 0 

Numerical solution 
computed from solving 
equation (21) 

10 


4.744 - 

0.268i 

4.772 - 0 

.244i 

11 


1.790 - 

3.354i 

1.615 - 3 

.226i 

12 


8.182 - 

.188i 

8.197 - 

.168i 

13 


11.450 - 

-141i 

11.461 - 

.126i 

14 


14.663 - 

.112i 

14.672 - 

.lOOi 

15 


17.851 - 

.093i 

17.858 - 

.082i 

16 


21.024 - 

.080i 

21.031 - 

.070i 

17 


24.190 - 

.070i 

24.196 - 

.060i 

18 


27.350 - 

.062i 

27.356 - 

.052i 

19 


30.505 - 

.055i 

30.512 - 

.045i 


TABLE VI.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 


(K = 5, B = 0.34 - 0.591) 


Radial 
number , 

mode 

£ 

Exact solution to 
equation: A^ 

transcendental 

- 5iB J2^^£^ " ° 

Numerical solution 
computed from solving 
equation (21) 

20 


2.157 - 

2.9301 

2.003 - 2 

.7701 

21 


6.212 - 

.2641 

6.234 - 

.2231 

22 


9.658 - 

.1691 

9.671 - 

.1521 

23 


12.940 - 

.1281 

12.949 - 

.1151 

24 


16.164 - 

.1041 

16.171 - 

.0921 

25 


19.360 - 

.0871 

19.366 - 

.0771 

26 


22.540 - 

.0751 

22.546 “ 

.0661 

27 


25.711 - 

.0661 

25.717 - 

.0571 

28 


28.875 - 

.0591 

28.881 - 

.0501 

29 


32.035 - 

.0531 

32.042 - 

.0431 
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TABLE VII.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE -WAVE SOURCE 


(K =10, B = 0.21 - 0.451) 


Radial 

number. 

mode 

£ 

Exact solution to 
equation: (X^) - 

transcendental 
lOlB ° 

Numerical solution 
computed from solving 
equation (21) 

0 


2.865 - 

0.194i 

2.912 - 0 

.185i 

1 


6.395 - 

.219i 

6.439 - 

.191i 

2 


9.736 - 

.179i 

9.769 - 

.152i 

3 


2.086 - 

5.024i 

1.742 - 4 

.713i 

4 


12.988 - 

.145i 

13.013 - 

.121i 

5 


16.198 - 

-121i 

16.220 - 

.lOOi 

6 


19.387 - 

.103i 

19.406 - 

.0841 

7 


22.573 - 

.090i 

22.580 - 

.0721 

8 


25.730 - 

.079i 

25.748 - 

.0621 

9 


28.892 - 

.071i 

28.911 - 

.0521 


TABLE VIII.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 

(K =10, 6 = 0.21 - 0.451) 


Radial 

number. 

mode 

£ 

Exact solution to 
equation: “ 

transcendental 
lOiB J^CX^) = 0 

Numerical solution 
computed from solving 
equation (21) 

10 


4.521 - 

0.2321 

4.571 - 0.2111 

11 


8.011 - 

.2021 

8.048 - .1751 

12 


2.142 - 

4.9271 

1.816 4.6181 

13 


11.321 - 

.1621 

11.349 - .1381 

14 


14.561 - 

.1331 

14.583 - .1121 

15 


17.766 - 

.1121 

17.785 - .0911 

16 


20.952 - 

.0961 

20.969 - .0791 

17 


24.126 - 

.0841 

24.142 - .0681 

18 


27.293 - 

.0751 

27.309 - .0591 

19 


30.455 - 

.0681 

30.473 ^ .0501 
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TABLE IX.- EIGENVALUES AT OPTIMUM ADMITTANCE FOR PLANE-WAVE SOURCE 


(K - 10, 3 = 0.21 - 0.451) 


Radial 

number. 

mode 

1 

Exact solution to 
equation: 

trans cendental 
- lOiB = 0 

Numerical solution 
computed from solving 
equation (21) 

20 


5.994 - 

0.2391 

6.041 - 0 

.2131 

21 


9.504 - 

.1881 

9.536 - 

.1621 

22 


2.323 - 

4.6281 

2.019 - 4 

.2911 

23 


12.822 - 

.1501 

12.846 - 

.1281 

24 


16.069 - 

.1241 

16.089 - 

.1041 

25 


19.280 - 

.1051 

19.297 - 

.0871 

26 


22.472 - 

.0911 

22.487 - 

.0751 

27 


25.651 - 

.0801 

25.666 - 

.0651 

28 


28.822 - 

.0721 

28.836 - 

.0571 

29 


31.987 - 

.0651 

32.003 - 

.0481 


TABLE X.- EIGENVALUES 


FOR TWO-STRIP PERIPHERALLY SEGMENTED LINER 


(T = 1, S = 2, K 


2, 3i = 1.54 - 1.291, 32 


0.0 + 0 . 01 ) 


M 


1 

3 


5 


1.363 

— 

1.9181 

1.346 

- 

0.3391 

1.276 

- 

0.2541 

3.450 

- 

.3881 

2.484 

- 

1.2831 

2.826 

- 

3.0441 

6.824 

- 

.2091 

2.830 

- 

3.0441 

3.323 

- 

.3741 

10.044 

- 

.1441 

3.416 

- 

.2661 

3.710 

- 

.4711 

13.225 

- 

.1091 

5.041 

- 

.2131 

3.638 

- 

2.3691 

16.391 

- 

.0881 

6.483 

- 

.2311 

4.616 

- 

.6841 

19.549 

- 

.0741 

6.813 

- 

.1881 

5.045 

- 

.2091 

22.703 

- 

.0631 

8.370 

- 

.1531 

6.469 

- 

.2301 

25.853 

- 

.0551 

9.831 

- 

.1501 

6.810 

- 

.1901 

29.002 

- 

.0481 

10.039 

- 

.1341 

7.810 

- 

.1801 


1.285 

2.823 

3.300 

3.581 

3.817 

5.044 

5.024 

5.765 

6.483 

6.813 


6 


10 


0.2301 

1.289 - 

0.2281 

3.0461 

2.825 - 

3.0481 

.3261 

3.295 - 

.3251 

2.3931 

3.599 - 

2.4141 

.3111 

3.830 - 

.3231 

.2121 

5.045 - 

.2121 

1.5521 

5.179 - 

1.5591 

.4091 

5.851 - 

.6281 

.1951 

6.514 - 

.1901 

.1831 

6.817 - 

.1791 
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TABLE XI.- EIGENVALUES FOR TWO-STRIP PERIPHERALLY SEGMENTED LINER 


(T = 1, S = 2, K = 5, = 1.54 - 1.29i, $2 = 0.0 + O.Oi) 


£ 







M 








1 



3 


5 

6 

10 

0 

2.815 

_ 

0.3881 

1.352 

_ 

0.1241 

1.319 - 0.1421 

1.258 

- 

0.1261 

1.281 

- 

0.1131 

1 

3.346 

- 

3.7861 

3.159 

- 

.1471 

3.111 - .1491 

3.083 

- 

.1481 

3.088 

- 

.1431 

2 

6.478 

- 

.4421 

4.786 

- 

.1581 

3.971 - .2041 

3.765 

- 

-2381 

3.869 

- 

.1931 

3 

9.828 

- 

.3141 

3.681 

- 

3.5791 

4.796 - .1511 

4.793 

- 

.1571 

4.796 

- 

.1531 

4 

13.067 

- 

.2441 

6.108 

- 

.4421 

4.561 - 3.0241 

5.404 

- 

.3551 

5.651 

- 

.3451 

5 

16.266 

- 

.1981 

6.516 

- 

.2151 

6.101 - .2901 

6.211 

- 

.2471 

6.330 

- 

.1861 

6 

19.446 

- 

.1661 

5.950 

- 

6.9381 

6.516 - .2111 

6.533 

- 

.2121 

6.566 

- 

.1931 

7 

22.614 

- 

.1411 

8.130 

- 

.1981 

5.959 - 6.9441 

5.960 

- 

6.9391 

5.959 

- 

6.9451 

8 

25.776 

- 

.1211 

9.611 

- 

.3151 

6.362 - 6.6241 

6.477 

- 

5.5631 

6.377 

- 

6.6341 

9 

28.934 

- 

.1031 

9.816 

- 

.2191 

7.544 - .2051 

6.392 

- 

6.5671 

7.557 

- 

.1821 


TABLE XII.- EFFECT INCREASING NUMBER OF STRIPS HAS ON EIGENVALUES OF 
PERIPHERALLY SEGMENTED LINER (T = 1 , K = 2, 

= 1.54 - 1.291, 62 = 0.0 + O.Oi) 


£ 


S 


2 

4 

10 

0 

1.289 - 0.2281 

2.046 - 0.7681 

1.568 - 1.7701 

1 

2.825 - 3.0481 

2.504 - .6781 

2.224 - 1.7151 

2 

3.295 - .3251 

2.982 - 2.8081 

2.978 - 1.5391 

3 

3.599 - 2.4141 

3.393 - 2.0131 

3.476 - .4791 

4 

3.830 - .3231 

3.912 ~ .2651 

3.917 - 1.3131 

5 

5.045 - .2121 

5.063 - .2751 

4.888 - 1.2081 

6 

5.179 - 1.5591 

5.098 - 1.2021 

5.138 - .3551 

7 

5.851 - .6281 

5.926 - 1.0331 

5.813 - .7041 

8 

6.514 - .1901 

6.529 - .1731 

6.524 - .2831 

9 

6.817 - .1791 

6.836 - .1671 

6.844 - .2521 
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TABLE XIII.- EFFECT INCREASING NUMBER OF STRIPS HAS ON EIGENVALUES OF 


PERIPHERALLY SEGMENTED LINER (T = 1, K = 5, 
3i = 1.54 - 1.291, ^2 = 0-0 + O.Oi) 




S 


£ 

2 

4 

10 

0 

1.281 - 0.113i 

1.869 - 0.332i 

2.603 - 0.453i 

2 

3.088 - .143i 

2.885 - .351i 

3.895 - .768i 

3 

3.869 - .193i 

4.154 - .3941 

4.882 - .663i 

4 

4.796 - .153i 

4.697 - .248i 

5.656 - 1.024i 

5 

5.651 - .345i 

6.077 - .862i 

4.732 - 3.420i 

6 

6.330 - .186i 

6.296 - .226i 

5.982 - .953i 

7 

6.566 - .193i 

6.664 - .302i 

5.625 - 3.1901 

8 

5.959 - 6.945i 

6.814 - l.Olli 

6.396 - .8091 

9 

6.377 - 6.634i 

6.092 - 6.825i 

6.864 - .6041 

10 

7.557 - .182i 

6.381 - 6.390i 

6.576 - 4.1651 


TABLE XIV.- REDISTRIBUTION OF ENERGY INTO MULTIPLE STANDING WAVES BY 


PERIPHERALLY SEGMENTED LINER (T = 1, K = 2, = 1 . 54 - 1.291, 

&2 = 0.0 + 0 . 01 ) 


Series 

coefficients 




Uniform liner 
values 


Peripheral liner 
with 2 strips 
(S = 2) 


Peripheral liner 
with 4 strips 
(S = 4) 


Peripheral liner 
with 10 strips 
(S = 10) 


^000 

^010 

^020 

^030 

^040 

^050 

^060 

^070 

^080 

^090 


0.213 

+ 

0.1881 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 

.000 

+ 

.0001 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 

.000 

+ 

.OOOi 


-0.347 

- 

0.1881 

.834 

+ 

.OOOi 

-.200 

+ 

.2081 

.187 

- 

.0281 

.001 

- 

.0841 

.067 

- 

.0841 

.014 

- 

.0271 

-.027 

+ 

.0311 

-.010 

- 

.OlOi 

.014 

- 

.0211 


-0,063 

+ 

0.2661 

-.340 

- 

.3361 

.539 

+ 

.OOOi 

-.275 

- 

.3471 

.077 

- 

.0751 

.155 

+ 

.1971 

-.006 

+ 

.0091 

.005 

- 

.0361 

-.018 

+ 

.0451 

-.018 

+ 

.0081 


0.525 

+ 

O.OOOi 

-.301 

+ 

.3371 

-.238 

- 

.0941 

-.015 

+ 

.1241 

-.115 

- 

.0241 

.004 

+ 

.0871 

-.086 

+ 

.0171 

.015 

+ 

.0941 

-.092 

- 

.0341 

.039 

+ 

.2721 






Nj|=f 


Figure 2.- Semi- infinite circular duct, 



Peripheral 

liner 


Figure 3.- Peripheral liner with periodicity of T = 4. 
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